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Abstract 

We review the closed time path formalism of Schwinger using a path integral approach. 
We apply this formalism to the study of pair production from strong external fields as 
well as the time evolution of a nonequilibrium chiral phase transition. 
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INTRODUCTION 



In 1961 in his classic paper "Brownian Motion of a Quantum Particle" 1 , Schwinger 
solved the formidable technical problem of how to use the action principle to study 
initial value problems. Previously, the action principle was formulated to study only 
transition matrix elements from an earlier time to a later time. The elegant solution 
of this problem was the invention of the closed time path (CTP) formalism. This 
formalism was first used to study field theory problems by Mahanthappa and Bakshi 1 . 

With the advent of supercomputers, it has now become possible to use this formal- 
ism to numerically solve important field theory questions which are presented as initial 
value problems. Two of these problems we shall review here. They are 

1. The time evolution of the quark-gluon plasma 2 . 

2. Dynamical evolution of a non-equilibrium chiral phase transition following a rel- 
ativists heavy ion collision 3 . 

The basic idea of the CTP formalism is to take a diagonal matrix element of the 
system at a given time t — and insert a complete set of states into this matrix element 
at a different (later) time t'. In this way one can express the original fixed time matrix 
element as a product of transition matrix elements from to t' and the time reversed 
(complex conjugate) matrix element from t' to 0. Since each term in this product is 
a transition matrix element of the usual or time reversed kind, standard path integral 
representations for each may be introduced. If the same external source operates in 
the forward evolution as the backward one, then the two matrix elements are precisely 
complex conjugates of each other, all dependence on the source drops out and nothing 
has been gained. However, if the forward time evolution takes place in the presence of 
one source J + but the reversed time evolution takes place in the presence of a different 
source J_, then the resulting functional is precisely the generating functional we seek. 
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Fig. 1 

Complex time contour C for the closed time path propagators. 
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so that, for example, 
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is a true expectation value in the given time-independent Heisenberg state \in). Here 
(p(t) = (p(x, t) and we are supressing the coordinate dependence and the integration 
over the spatial volume in what follows for notational simplicity. 

Since the time ordering in eq. (0) is forward (denoted by T) along the time path 
from to t' in the second transition matrix element, but backward (denoted by T*) 
along the path from t' to in the first matrix element, thus the name: closed time 
path generating functional. If we deform the backward and forward directed segments 
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of the path slightly in opposite directions in the complex t plane, the symbol Tc may 
be introduced for path ordering along the full closed time contour, C depicted in Fig.l. 

This deformation of the path corresponds precisely to opposite ie prescriptions along 
the forward and backward directed segments, which we shall denote by C± respectively 
in the following. 

If we have an arbitrary initial density matrix p then we have instead: 



Z[J+,J_,p] = Trip T*exp 



t' 

-i i dtJ-(t)<p(t) 
o 



Texp 



a 

i i dtJ+(t)<p(t) 
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[XV][2V][Xty] (vW) ( V '\T*ex P 
t' 



f 



-% I dtJ_(t)(f)(t) 
o 



x (-01 Texp 



i l dtJ + (t)(f)(t) 
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\<P) 



Variations of this generating function will yield Green's functions in the state specified 
by the initial density matrix, i.e. expressions of the form, 

Tr{p<^i)0(*2)</>(i 3 )---}- (4) 

Introducing the path integral representation for each transition matrix element in 
eq. (^|) results in the expression, 

z[J+,J-,p] = jpvW] (^W) J[v^f\v<t> + ]jyv^\ x 

exp 



oo 

il dt(L[ 
o 



-L\ 



J4 



J_ 



(5) 



where L is the classical Lagrangian functional, and we have taken the arbitrary future 
time at which the time path closes t' — > 00. 

The double path integral over the fields + and 0_ in (|) suggests that we introduce 
a two component contravariant vector of field variables by 



with a corresponding two component source vector 



r 



1,2 



a = 1,2 . 



(6) 



(7) 



Because of the minus signs in the exponent of (f|), it is necessary to raise and lower 
indices in this vector space with a 2 x 2 matrix with indefinite signature, namely 



(3) 



c ab = diag (+1, -1) = c' 



ab 
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so that, for example 



J a c nh $ h 



j ab w = J + (j) + - J_0_ . (9) 
These definitions imply that the correlation functions of the theory will exhibit a matrix 



structure in the 2x2 space. For instance, the matrix of connected two point functions 
in the CTP space is 

G ab {t,t') 



5 2 W 



5J a (t)5J b (t') 

Explicitly, the components of this 2x2 matrix are 



j=o 



(10) 



G 21 (M') 
G 12 (M') 
G u (t,t') 

G 22 (t,t f ) 
Notice that 



G>(M')= *Tr{p <5>(t)$(t')} con , 
G<(t,t') =iTr{ P WMt)}con 

iTr\p T[$(i)$(t')]} = Oit^G^t') + Oit'^G^t.t 1 ) 

t i con 

iTr {p T* [$(*)$(*')]} = d{t', t)G>(t, + 0(t, t , )C<(t, t') 

>. J con 



G n (i,t) =G™(t,t) 



22, 



with the usual convention that 



6(t,t) 



121 



(13) 



The 2x2 matrix notation originated with Schwinger's classic article in I960 1 . 

In what follows we will use an alternative generating functional 4 using the Complex 
Path Ordered Form: 



[ViJ)]{(p'\T* exp 



dtJ-(t)<j>(t) 



\ip) (^|Texp 



dt J + (t)(p(t) 



\<p) 



((^'iT^exp i J dt J{t)4>(t) 



\<f) 



(14) 



so that @ may be rewritten more concisely in the CTP complex path ordered form, 



Zc[J,p) 



Tr p T c exp 



% \ dtJ(t)(/>(t) 



pV] J [Dip 2 ] (ip l \p\ip 2 ) J\ [X?0]exp Ujdt (L[<f>] + J0) 



• (15) 



This is identical in structure to the usual expression for the generating functional 
in the more familiar in-out formalism, Only difference - path ordering according to the 
complex time contour C replacing the ordinary time ordering prescription along only 
C + . For example, the propagator function becomes 



G(t,t') = # c (M0G>(M0 + Wt)G<(M') 
= 6c(t,t')G 21 (t,t>) + d c (t',t)G 12 (t,t') 



(16) 
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where 8c is the CTP complex contour ordered theta function defined by 



9 e (t,t') 



8(t,t') for t, t' both on C + 

9(t', t) for t, t! both on C_ 

1 for t on C_ , t' on C + 

for i on C + , i' on C_ 



With this definition of G(t,t') on the closed time contour, the Feynman rules are the 
ordinary ones, and matrix indices are not required. In integrating over the second half 
of the contour C_ we have only to remember to multiply by an overall negative sign to 
take account of the opposite direction of integration, according to the rule, 



dt = dt- dt . (18) 
c Jo C+ JO C- 

A second simplification is possible in the form of the generating functional of (|I5D, 
if we recognize that it is always possible to express the matrix elements of the density 
matrix as an exponential of a polynomial in the fields 5 

(<p l \p\<p 2 ) =exp[R + R a (to)cp a (to) + Rab(to)<p a (to)<p b (t ) + ...} . (19) 

Since any density matrix can be expressed in this form, there is no loss of generality 
involved in expressing p as an exponential. If we add this exponent to that of the 



action in (15[), and integrate over the two endpoints of the closed time path ip 1 and ip 2 , 
then the only effect of the non-trivial density matrix p is to introduce source terms into 
the path integral for Z C [J, p] with support only at the endpoints. This means that the 
density matrix can only influence the boundary conditions on the path integral at t — 0, 
where the various coefficient functions R a , R a b, etc. have the simple interpretations of 
initial conditions on the one-point (mean field), two-point (propagator), functions etc. 
It is clear that the equations of motion for t ^ are not influenced by the presence 
of these terms at to = 0. In the special case that the initial density matrix describes 
a thermal state, pp = exp{— f3H} then the trace over pp may be represented as an 
additional functional integration over fields along the purely imaginary contour from 
t = —i/3 to t = traversed before C_ in Fig. 1. In this way the Feynman rules for real 
time thermal Green's functions are obtained 6 . Since we consider general nonequilibrium 
initial conditions here we have only the general expression for the initial p above and 
no contour along the negative imaginary axis in Fig. 1. 

To summarize, we may take over all the results of the usual scattering theory gen- 
erating functionals, effective actions, and equations of motion provided only that we 

1. substitute the CTP path ordered Green's function(s) (|I6|) for the ordinary Feyn- 
man propagators in internal lines; 
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2. integrate over the full closed time contour, C, according to fll8|) ; and 

3. satisfy the conditions at t — corresponding to the initial density matrix p. 



Closed time path contour and causality 

Rules for evaluating the time integrals using the closed time path [CTP] contour 
shown in Fig. (1). 

The integration path is given by 

r- roo roc 

dt = dt - / dt. (20) 

Jc J0-.C+ JQ:c- 

The causal Green's functions are given by functions of the form, 

A(t,f) = e c (t,f)^>(t,f) + e^t'^A^tX) , (21) 

where O c (t,t') is defined in eq. (|lTD These causal Green's functions are symmetric. 

A > {t,t')=A < {t',t) 

To prove causality of any graph we need two lemmas. 

1. Lemma 1- a loop of two causal functions, such as self energy graph, is another 
causal function. 

B(t,t') = e^MO^MO + ec^t^M') 

c(m') = e^o^^O + e^o^^O, (22) 

the self energy graph 

A(t,t') = tB(t,t')C(t,t'), 

then 

A(t, = 9 c (t, ti)A> (t, t') + 9 c (t', t)A < (t, t') (23) 

where 

A > , < (t,t')=tB >t< (t,t')C > , < (t,t') (24) 

2. Lemma 2 -Matrix product of two causal functions is causal 

A(h,t 3 ) = [ dt 2 B(t!,t z )C(t 2 , t 3 ) , (25) 

J c 
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we find then 

A(t,t') = e c (M'M>(M') + e c (t , ,t)A < (t,0 , (26) 

where 

A>(ti,* 3 ) = - F dt 2 B > (t 1 ,t 2 ) [C>(f 2 ,t 3 ) -C<(t 2 ,t 3 )] 
< Jo < 

+ / 1 d* 2 [B > (ti,* 2 )- J B < (*i,t 2 )]C>(* 2 ,* 3 ). (27) 
jo < 

Now consider the product of three causal functions: 

A(t u U) = j &t 2 j &hB{t u t 2 )C{t 2 ,h)D{hM) • (28) 

We can work this case out by applying the second lemma from left to right. That is, 
we can let 

E(t u t 3 ) = jdt 2 B(t u t 2 )C(t 2 ,t 3 ) . (29) 
Then E(t 1: t 3 ) is causal. We are then left with: 

A(t l} t A ) = JdtsEfat^DfaU) , (30) 

and so A is also causal. 

After doing the integrals sequentially one is left with 

/(*) = jdhF{tM) = /V>(Mi) - ^<(Mi)], (31) 

which explicitly displays the causality (dependence only on earlier times). 

To see how these rules work in practice, consider the terms contributing to order 
1/N to the induced current determining the backreaction on an initially strong electric 
field 4 . The current is just tr^^G}. Where G is the full fermion propagator to order 
1/N. 

Using the above lemmas we obtain that the Maxwell eqs. of motion take the form, 

d„F^(x) = (f(x)) = -—tT{T[G > (x,x)+G < (x,x)]} 

2e 2 /■* /■*! f 

+ — / dhd 3 ^ dt 2 d 3 x 2 lmtr\iv \G > (x,x 1 ) - G < (a;,a; 1 )] x 
N Jo jo I 

[E < (x 1 ,x 2 )G > (x 2 ,x) - E > (x 1 ,x 2 )G < (x 2 ,x)] J . (32) 

Here: 

T, < (x 1 ,x 2 ) = i^G < (x 1 ,x 2 )YD u ^ 1< (x 1 ,x 2 ). (33) 
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This immediately displays the causality of the result: that is only previous times 
contribute to the space-time integrations. 

TIME EVOLUTION OF THE QUARK GLUON PLASMA 

Our model for the production of the quark-gluon plasma begins with the creation 
of a flux tube containing a strong color electric field. We assume the kinematics of 
ultrarelativistic high energy collisions results in a boost invariant dynamics in the lon- 
gitudinal (z) direction where all expectation values are functions of the proper time 
t = \Jt 2 — z 2 .We therefore introduce the light cone variables r and rj, which will be 
identified later with fluid proper time and rapidity . These coordinates are defined 
in terms of the ordinary lab-frame Minkowski time t and coordinate along the beam 
direction z by 

^ = rsinh?7 , t = rcosh?7. (34) 
The Minkowski line element in these coordinates has the form 

ds 2 = -dr 2 + dx 2 + dy 2 + T 2 dif . (35) 

Hence the metric tensor is given by 

^ = diag(-l,l,l,r 2 ). (36) 

For simplicity, here we discuss pair production from an abelian Electric Field and the 
subsequent quantum back-reaction on the Electric Field. The physics of the problem 
can be understood for constant electric fields as a simple tunneling process. If the 
electric field can produce work of at least twice the rest mass of the pair in one compton 
wavelength, then the vacuum is unstable to tunnelling. This condition is: 

— > 2mc 2 . (37) 
mc 

The problem of pair production from a constant Electric field (ignoring the back 
reaction) was studied by J. Schwinger in 1951 7 . The WKB argument is as follows: 
One imagines an electron bound by a potential well of order \Vq\ ~ 2m and submitted 
to an additional electric potential eEx . The ionization probability is proportional to 
the WKB barrier penetration factor: 

exp[-2 [ V ° /e dx{2m(V - \eE\x)} 1/2 } = exp(-^m 2 /\eE\) (38) 
Jo 3 
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In his classic paper Schwinger was able to analytically solve for the effective Action in 
a constant background electric field and determine an exact pair production rate: 



w = \a 



£7(2vr 2 )]£- 



\n+l 



n=l 



exp ( — rnim 2 / \ eE | ) . 



(39) 



By assuming this rate could be used when the Electric field was slowly varying in 
time, the first back reaction calculations were attempted using semi classical transport 
methods. Here we use the CTP formalism and perform the field theory calculation. 
The lagrangian density for QED in curvilinear coordinates gives rise to the action 



S = J d d+l x (det V) ^#7*%* + |(vj»f^ - imff - -F^F^ 



where 8 



(40) 
(41) 



= (d^ + T^-ieA^ 
From the action ([10) we obtain the Heisenberg field equation for the fermions, 

(7 M V M + m) * = , (42) 

which takes the form 



* = 



(43) 



Variation of the action with respect to A u yields the Maxwell equations: If the electric 
field is in the z direction and a function of r only, we find that the only nontrivial 
Maxwell equation is 



1 dE(r) 
t dr 



2r 



(44) 



We expand the fermion field in terms of Fourier modes at fixed proper time r, 

= J [dk] ^[6 s (k)^ s (r)e^e i P- x + ^-k^Z^r)^^-*-*]. (45) 



The ijj^ s then obey 



7 



dr 2t 



V£(r) = 0, 



We square the Dirac equation by introducing 



ks 



27 ir v + m 



Xs 



ks 



(46) 



(47) 
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The spinors Xs are chosen to be eigenspinors of 7°7 3 , 

7°7 3 Xs = KXs 

with A., = 1 for s = 1, 2 and X s = — 1 for s = 3,4. They are normalized, 

Xr Xs 25 rs . 



(4f 



(49) 



The sets s = 1,2 and s = 3,4 are two different complete sets of linearly independent 
solutions of the Dirac equation Inserting (47]) into the Dirac equation fl46|) we obtain 
the quadratic mode equation 



(^+<4-<A.*,)y£(r) = 0, 



where now 
We obtain 



1 dE(r) 
r c?r 



2e 

^2 



E /[rfk](ki + m 2 )A s |/, 

s=l J 



ksl ' 



(50) 



(51) 



(52) 



We would like to compare the field theory calculation with a semiclassical transport 
approach which uses as a source of particle production Schwinger's production rate for 
time independent fields. 

Assuming boost-invariant initial conditions for /, invariance of the Boltzmann- 
Vlasov assures that the distribution function is a function only of the boost invariant 
variables (r, 77 — y) or (r,p v ). The kinetic equation reduces to 



df „ , v df 



±[l±2/(p,r)]er|S(r) 



x In 



1 ± exp 



7r(m 2 + p 2 ± ) 
e\E(r)\ 



Turning now to the Maxwell equation, we have that 

_t- — A . — Acond 1 -pol 

dr h ^ +3r < ' 



(53) 



(54) 



where j cond is the conduction current and j v ° is the polarization current due to pair 
creation 9 . 
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Fig. 2 

Proper time evolution of the electric field E(t) and the fermion current J^(t) for 
initial conditionsi?(Y = 1) = 4.0. The field theory calculation is compared to the 
semiclassical transport approach with a Schwinger source term. 



Thus in (1 + 1) dimensions we have 
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■pol 
'V 



2_ r dp, Df 



F Tr i J 2-ktp/ Dt 



±[l±2/(p„ = 0,r)] 



mer 



sign[_E(r)] In 1 ± exp 




7T 



\eE(r)\)_ 



(55) 



In figure two we compare the results of a numerically solving the back reaction equa- 
tions in the field theory with the semiclassical transport approach. We find that for this 
approximation (no rescattering of quarks), the semiclassical method does reasonable 
well when compared to a coarse grained in time and momentum field theory calculation. 

DYNAMICAL EVOLUTION OF A NON-EQUILIBRIUM CHIRAL PHASE 
TRANSITION 

Recently there has been speculation that disoriented chiral condensates (DCC's) 
can be formed following a heavy ion collision, and these condensates, formed during a 
quenched phase transition from the unbroken high temperature phase, could lead to 
events with a nonequilibrium distribution of charged to neutral pions. To see whether 
these ideas made sense we studied numerically the time evolution of pions produced 
following a heavy ion collision using the linear sigma model, starting from the unbroken 
phase. The quenching (if present) in this model is due to the expansion of the initial 
Lorentz contracted energy density by free expansion into vacuum. Starting from an ap- 
proximate equilibrium configuration at an initial proper time r in the disordered phase 
we study the transition to the ordered broken symmetry phase as the system expands 
and cools. We determined 3 the proper time evolution of the effective pion mass, the 
order parameter < a > as well as the pion two point correlation function. We studied 
the phase space of initial conditions that lead to instabilities (exponentially growing 
long wave length modes) which can lead to disoriented chiral condensates. The model 
we used to study the chiral phase transition is the a model in the large N approxima- 
tion. This model has the correct chiral properties and gives a reasonable description of 
low energy pion dynamics. The 0(4) a model is described by the Lagrangian: 




<9$ - -A($ • $ - v 2 ) 2 + Ha}. 



(56) 



The mesons form an 0(4) vector 



$ = (a, 7Ti) 



Introducing the order parameter: 



X = A($ • $ - v 2 ). 
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we have the alternative Lagrangian: 

U = + + ^ + \xv 2 + Ha (57) 

Perform the Gaussian path integral over the $ field. Evaluate the remaining x 
integral at the stationary phase point . Legendre transforming: 

r[$, X ] = J d 4 x[L 2 ^, X ,H) + ^iVtr InGo 1 ] 

G^(x,y) = z[n + X (x)} 5\x-y) 
p + x(x)]Tn = [n + x (x)]a = H 



x = -\v 2 + A(cr 2 + 7T ■ 7r) + XNG (x,x). (58) 

As in the quark-gluon plasma problem we consider the kinematics of an ultrarelativistic 
Heavy Ion Collision which possesses longitudinal Boost invariance as the center of mass 
energy goes to infinity. Energy densities become function of the proper time only. 
Natural coordinates are the proper time r and the spatial rapidity rj defined as 

r^(t 2 -x 2 )^, ^ilog(^). (59) 

We assume that the mean (expectation) values of the fields $ and x are functions of r 
only (homogeneity in the constant r hypersurface) 



r- 1 d T rd T $,(r) + X (r) $,(r) = H6 a 

x (t) = A(-t; 2 + $ J 2 (r) + iV<0 2 ( a ;,r)>). (60) 

On the otherhand the fluctuation fields (which are quantum operators) are functions 
of both x and t and obey the sourceless equation: 

{r^drrdr -r- 2 d 2 v -dl + x(x))<P(x,T)=0. (61) 



G (x, y; r) =< <f>(x, r) <j>{y, r) > . (62) 
We expand this field in an orthonormal basis: 

4>(i],x ± ,t) = J[d 3 k](exp(ikx)f k (r) a k + h.c.) (63) 
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where kx = k v i] + k±xj_, [d 3 k] = dk v d 2 k±/(2'ir) 3 and the mode functions /k(T) evolve 
according to (a dot here denotes the derivative with respect to the proper time r): 

A+(^ + kl + x(r) + ^- 2 )h = 0. (64) 

X (r) = \(-v 2 + **(r) + M 3 k]|/ k (r)| 2 (1 + 2 n k )). (65) 

If we assume the initial density matrix is one of local thermal equilibrium then we 
have at r = r (the surface of constant energy density and temperature T ) that: 

1 

Uk ~ e PoEfi k Z i 

where f3 Q = 1/T and = \Jk 2 + x(r ). 

In choosing initial conditions we assumed that the initial value of \ was determined 
by the equilibrium gap equation for an initial temperature of 200MeV. The phase 
transition in this model occurs at a critical temperature of 160MeV. We chose initial 
u and 7T, expectation values consistent with the constraint 

7f 2 (r ) + a 2 (r ) = 4 (66) 

where ot is the equilibrium value of $ at the initial temperature T. which we choose to 
be a temperature of 200MeV '. We varied the value of the initial proper time derivative 
of the sigma field expectation value and found that there is a narrow range of initial 
values that lead to the growth of instabilities. Namely 

.25 < \a\ < 1.3 (67) 

Surprisingly when \&\ > 1.3 instabilities no longer occur. 

Figures 3-4 summarize the results of the numerical simulation for the evolution of 
the system (|65|)-(|64]). We display the auxiliary field x i n units of fm~ 2 , the classical 
fields $ in units of /m _1 and the proper time in units of /m _1 (IfrrT 1 = 197 MeV). 
In Fig. 3 the proper time evolutions of the auxiliary field \ field is presented for four 
different initial conditions. 

We are interested in knowing how our results differ from the case where the system 
evolves in local thermal equilibrium which is described by two correlation lengths, the 
inverse of the effective pion mass associated with x, and the inverse of the proper time 
evolving effective temperature T(r) = T (^)3 discussed earlier. We see from Fig. 4 
that in the case that cr(l) = <J T , 7r*(l) = and <j(l) = —1, where maximum instability 
exists, complex structures are formed as contrasted to the local thermal equilibrium 
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evolution. The interpolating phase space distribution n(k v , kj_, r) obtained numeri- 
cally, clearly exhibits a larger correlation length in the transverse direction than the 
equilibrium one and has correlation in rapidity of the order of 1-2 units of rapidity. We 
notice that in both directions there is structure which does not lend itself to a simple 
interpretation. On the other hand the local thermal equilibrium evolution is quite reg- 
ular apart from the normalization of the distributions that are changing with time due 
to oscillation in the quantity x( r ) which is damped to its equilibrium value once the 
system expands sufficiently. 
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Proper time evolution of the x fi e ld for four different initial conditions with f n 
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